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I .  INTRODUCTION 


Experimental  values  of  the  attenuation  of  acoustic  waves  below 

100  Hz  have  steadily  accumulated  during  the  last  15  years.  Kibblewhite 

and  Hampton1  have  summarized  and  classified  the  available  data  and 

provided  a  bibliography.  The  present  report  is  concerned  with  the 

excess  attenuation  values  below  100-200  Hz.  Kibblewhite  and  Hampton 

introduce  a  term  a  called  the  scattering  coefficient,  expressed  in 
s 

dB/km,  as  a  measure  of  this  excess  attenuation.  Their  Table  I  suggests 
that  this  coefficient  is  related  to  the  water  mass  in  which  the  acoustic 
profile  is  located. 

It  is  not  possible  at  the  present  time  to  propose  a  mechanism  for 
this  scattering.  Nonetheless,  it  seems  profitable  to  assign  various 
depth  and  frequency  characteristics  to  this  scattering  mechanism  and  to 
make  quantitative  predictions  of  attenuation  for  long  distance  propaga¬ 
tion.  The  results  of  these  calculations  will  provide  significant  guid¬ 
ance  in  the  search  for  and  identification  of  the  relevant  scattering 
mechanisms.  These  calculations  extend  a  line  of  calculations  initiated 

p 

by  Guthrie  in  his  Ph.D,  dissertation.  Guthrie's  work  is  summarized 
briefly  by  Kibblewhite,  Shooter,  and  Watkins^  in  their  work  on  the 
attenuation  of  deep  water  ambient  noise. 
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II.  MODELING 


This  investigation  was  guided  by  our  experience  in  computing 
propagation  loss  in  shallow  water  propagation.  At  low  frequencies  sig¬ 
nificant  propagation  loss  results  from  the  acoustic  energy  that  propa¬ 
gates  in  the  bottom.  Normally  this  loss  would  be  expected  to  depend 
linearly  on  the  frequency. ^  However,  the  attenuation  of  clastic  sedi- 
ments  varies  significantly  with  the  nature  of  these  sediments.'  Since 
the  lower  frequency  waves  penetrate  more  deeply,  the  dependence  of  prop¬ 
agation  loss  on  frequency  may  be  more  strongly  influenced  by  the 
sedimentary  column  than  by  the  frequency  dependence  itself.  In  any  case 
these  considerations  suggest  that  one  should  not  rule  out  possible  depth 
or  frequency  distributions  on  a  priori  arguments  that  they  are 
implausible.' 

I 

The  properties  of  the  scatterers  in  the  deep  ocean  can  vary  with 
both  depth  in  the  water  and  frequency.  Their  effect  on  the  acoustic 
wave  will  depend  on  the  variation  of  the  acoustic  pressure  with  depth 
and  thence  on  the  velocity  structure  of  the  water  and  on  the  mode  number 
of  the  wave.  The  following  sequence  was  used  in  the  calculations. 
After  the  sound  velocity  profile  was  chosen,  the  eigenfunctions  for  this 
velocity  structure  in  the  absence  of  attenuation  were  oomputed  using  the 
NEMESIS  program. ^  A  specific  depth  distribution  of  attenuation  was 
assumdd  both  in  the  water  and  in  the  sediments  and  the  resulting  propa¬ 
gation  loss  was  computed  by  a  perturbation  technique  similar  to  that 
2 

used  by  Guthrie. 

Urick®  has  an  excellent  summary  of  attenuation  in  sea  water  and  of 
the  dependence  on  ooeanographio  parameters.  Figure  1  is  a  plot  of 

Q 

attenuation  versus  depth  based  on  experimental  data  of  Bezdek*  and 

oorrected  for  salinity  and  temperature  variation  using  formulas . by. 

Sohulkin  and  Marsh.10  This  ourve  shows  a  maximum  attenuation  near  the 
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sound  channel  axis  followed  by  a  linear  decrease  with  depth.  These 
curves  contrast  with  the  models  used  in  the  present  analysis  where  the 
attenuation  is  assumed  to  increase  all  the  way  to  the  surface.  This 
difference  is  not  significant  in  the  near-surface  layers  since  the 
acoustic  pressure  vanishes  at  the  surface  of  the  water. 

In  the  models  analyzed  below,  it  is  assumed  that  the  scatterers  are 
concentrated  in  the  upper  500  m  of  the  water  column.  This  decision  was 
influenced  in  part  by  oceanographic  data  such  as  the  illustrative 
examples  shown  in  Fig.  2.  Mellen  and  Browning11  show  that  pH,  and 
consequently  attenuation,  in  the  ocean  increases  strongly  near  the  sea 
surface.  Another  consideration  was  the  amplitude  versus  depth  curves 
for  first  order  acoustic  modes  in  the  deep  ocean.  Illustrative  curves 
for  the  amplitude  of  the  first  mode  for  12.5  and  100  Hz  are  shown  in 
Fig.  3.  These  curves  illustrate  the  concentration  of  the  acoustic 
energy  in  the  upper  1000  m.  They  also  show  how  the  12.5  Hz  mode  is  more 
strongly  affected  than  the  100  Hz  mode  when  the  attenuation  mechanisms 
are  concentrated  near  the  surface  of  the  ocean. 

Figures  4  and  5,  which  are  reproduced  from  Kibblewhite  and 
Hampton,1  show  the  phenomena  under  study.  In  Fig.  4  the  attenuation  has 
a  minimum  near  60  Hz  and  appears  to  increase  at  lower  frequencies. 
Figure  5  is  of  interest  in  that  it  shows  the  highest  observed  values  of 
attenuation  at  these  low  frequencies,  a  result  perhaps  of  the  arctic 
convergence  zone. 

Figures  6  and  7  provide  a  visual  comparison  of  the  observed 

attenuations  in  water  and  sediments  with  the  boric  acid  attenuation  of 

12 

Thorp  for  frequencies  of  12.5  and  100  Hz,  respectively.  It  is  note¬ 
worthy  that  at  12.5  Hz  the  attenuation  in  the  water  column  approaches 
that  of  sediments  and  is  greatly  in  excess  of  the  value  for  boric  acid. 
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FIGURE  3 

ACOUSTIC  PRESSURE  FOR  THE  FIRST  NORMAL  MODE  AT  12.5  Hz  (SOLID  CURVE) 

AND  100  Hz  (DASHED  CURVE) 
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FIGURE  4 

ATTENUATION  DATA  FOR  A  MID-LATITUDE  SOUND  VELOCITY  PROFILE  "TYPE  B" 
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FIGURE  5 
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ATTENUATION  DATA  FOR  THE  ARCTIC  CONVERGENCE  IN  THE  EASTERN  NORTH  PACIFIC  OCEAN 
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FIGURE  6 

COMPARISON  OF  OBSERVED  ATTENUATION  AT  12.5  Hz  IN  WATER  AND  SEDIMENTS 
AND  THE  VALUE  PREDICTED  BY  THE  THORP  FORMULA 
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FIGURE  7 

COMPARISON  OF  OBSERVED  ATTENUATION  AT  100  Hz  IN  WATER  AND  SEDIMENTS 
AND  THE  VALUE  PREDICTED  BY  THE  THORP  FORMULA 
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III.  NUMERICAL  PREDICTIONS  OF  SOUND  ATTENUATION 


The  arguments  outlined  above  suggest  a  provisional  use  of  the  three 
attenuation  profiles  shown  in  Fig.  8.  The  total  attenuation  coefficient 
o(z,f)  is  the  sum  of  the  value  <i^(f)  predicted  by  the  Thorp  curves  and 
the  contribution  of  a  frequency  dependent  profile  such  as  one  of  those 
shown  in  Fig.  8;  that  is,  the  attenuation  will  be  computed  from  the 
formula 

<>(z,f)  =  <*T(f)  +  k(z)f2  .  (1) 

Figure  9  shows  the  velocity  profile  used  in  the  propagation 
calculations.  The  water  depth  is  4883  m.  This  is  underlaid  by  50  m  of 
clay.  The  velocity  contrast  at  the  ocean  bottom  is  unity,  and  the  gra¬ 
dient  of  the  sound  speed  in  the  sediment  is  1.23  sec-1. 

As  mentioned  before,  the  normal  modes  are  computed  first  without 

14  5 

attenuation,  using  the  NEMESIS  program.  *  The  effect  of  attenuation, 
Eq.  (1),  is  introduced  oy  a  perturbation  technique.  These  complex 
normal  modes  are  then  used  to  calculate  propagation  loss  PL(r)  at  large 
ranges  r.  The  resulting  values  of  propagation  loss  are  converted  to 
attenuation  coefficients  with  the  formula 

n  =  (1/2,000)[PL(1 ,000)  -  PL(3,000)  -10  log( 3 ,000/1 ,000 )]  ,  (2) 

which  corrects  for  a  cylindrical  spreading  loss.  This  method  provides  a 
close  parallel  to  the  way  experimental  data  are  treated. 

Figure  10  shows  the  results  of  this  calculation  for  profiles  I 
and  II  of  Fig.  8.  In  each  of  these  curves  the  surface  value  of  attenua¬ 
tion  is  1  dB/km  at  100  Hz.  The  two  curves  labeled  T  and  0.5T  are  the 

Thorp  curves  for  two  different  concentrations  of  boric  acid.  The 
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FIGURE  9 

PROFILE  USED  IN  CALCULATING  ATTENUATION 
versus  FREQUENCY 
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FIGURE  10 

CALCULATED  COEFFICIENT  FOR  PROFILES  I  AND  II  WITH  A 
FREQUENCY  SQUARED  DEPENDENCE 

ARL:UT 
AS-81-657 
KCF - GA 
5 • 22  •  81 


THORP  CURVES  PRESENTED  AS  A  REFERENCE 


attenuation  coefficients  predicted  by  these  two  models  and  shown  in 
Fig.  10  are  not  acceptable.  Both  profiles  yield  excess  attenuation  but 
only  profile  II  shows  any  flattening  at  the  lower  frequencies. 

The  two  curves  shown  in  Fig.  10  are  calculated  for  source  and 
receiver  at  a  depth  of  600  m,  the  depth  of  the  sound  channel  axis. 
Since  some  data  are  measured  with  the  receiver  at  a  shallower  depth,  the 
attenuation  coefficient  for  profile  II  is  repeated  with  receiver  depth 
at  300  m.  Figure  11  shows  the  results  for  a  source  at  600  m  and 
receivers  at  300  and  600  ro  and  attenuation  profile  II.  As  might  be 
expected  the  predicted  attenuation  is  larger  at  300  m  but  the  frequency 
dependence  has  increased.  Further  calculations  with  frequency  squared 
attenuation  values  were  made  but  they  did  not  lead  to  the  flattening  of 
attenuation  values  observed  in  the  experimental  data. 

The  next  set  of  calculations  are  based  on  a  formula  similar  to 
Eq.  (1)  except  that  the  frequency  factor  f  in  the  last  term  is  replaced 
by  f  to  the  first  power.  These  calculations  are  made  with  the  attenua¬ 
tion  profile  II  of  Fig.  8  with  three  different  values  of  attenuation  at 
z=0.  Curves  labeled  IV,  V,  and  VI  correspond  to  a  surface  attenuation 
at  100  Hz  of  0.1  dB/km,  0.01  dB/km,  and  1.0  dB/km,  respectively.  The 
resulting  attenuation  coefficients,  shown  in  Fig.  12,  yield  a  more 
nearly  satisfactory  frequency  dependence.  Curve  VI  actually  shows  a 
minimum  as  do  some  of  the  experimental  curves,  but  the  magnitudes  of  the 
theoretical  curves  are  low.  It  is  interesting  to  note  that  curve  VI  of 
Fig.  12  agrees  very  closely  with  the  experimental  data  for  a  "tropical" 
sound  velocity  shown  in  Fig.  7  of  Kibblewhite  and  Hampton.1  Figure  13 
shows  a  second  calculation  of  case  VI  with  the  receiver  located  at  300  m 
in  addition  to  the  curve  for  600  m.  The  presence  of  the  minimum  is  seen 
to  be  receiver  depth  dependent. 

As  an  alternative  to  the  attenuation  profiles  shown  in  Fig.  8,  the 
formula 

rt(z,f)  =  oQ(f)  exp(-z/h)  +oT(f)  ,  (3) 
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FIGURE  12 

CALCULATED  ATTENUATION  COEFFICIENTS  FOR  PROFILE  II 
WITH  A  LINEAR  FREQUENCY  DEPENDENCE 
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is  tried  with  h  =  200  m  and  aQ(f)  =  0.1(f/100)  dB/km  where  f  is  in 
hertz.  The  two  curves  shown  in  Fig.  14,  computed  for  receiver  depths  of 
300  and  600  m,  do  not  show  the  desired  frequency  dependence. 

The  attenuation  coefficients,  selected  to  describe  a  profile,  need 
to  be  consistent  with  the  measured  values  at  frequencies  above  1  kHz  and 
in  particular  above  10  kHz.  Based  on  the  work  of  Schulkin  and  Marsh, ^ 
the  surface  attenuation  value  at  10  kHz  should  be  approximately  1  dB/km. 
The  surface  attenuation  value  for  a  linearly  frequency  dependent  term, 
necessary  in  matching  the  low  frequency  measurements,  would  exceed  this 
10  kHz  value  by  at  least  an  order  of  magnitude. 

The  behavior  of  the  attenuation  coefficients  computed  so  far  forces 

one  to  the  conclusion  that  a  successful  scattering  model  will  be  very 

insensitive  to  frequency.  Consequently  the  attenuation  function  of 

Eq.  (3)  is  tried  for  the  special  case  a  ( f)=constant.  This  is  a  model 

0  2 

of  the  attenuation  that  was  proposed  and  investigated  by  Guthrie.  The 

three  attenuation  curves  plotted  in  Fig.  15  are  computed  for 

a  =0.11  dB/km  and  for  h  =  200  m,  250  m,  and  500  m.  In  each  case  the 
o 

attenuation  has  a  minimum.  The  frequency  of  the  minimum  increases  as 
the  characteristic  depth  h  increases.  These  curves  reproduce  rather 
well  the  behavior  of  the  experimental  data  summarized  by  Kibblewhite  and 
Hampton1  and  the  two  examples  reproduced  in  Figs.  4  and  5  of  the  present 
report.  The  magnitude  of  the  attenuation  in  the  flat  portions  of  the 
curves  spans  very  nicely  the  values  of  ag  listed  in  Kibblewhite  and 
Hampton's  Table  I. 
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FIGURE  15 

CALCULATED  ATTENUATION  COEFFICIENTS  FOR  AN  ATTENUATION  PROFILE  WITH  AN 
EXPONENTIAL  DEPTH  DEPENDENCE  AND  INDEPENDENT  OF  FREQUENCY  , 

THE  CHARACTERISTIC  DEPTH,  H,  AS  A  PARAMETER  £L81 
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IV .  COMMENTS 


The  calculations  presented  here  do  not  furnish  a  mechanism  for  the 
attenuation  of  deep  ocean  acoustic  waves,  but  they  do  provide  a  set  of 
quantitative  characteristics  that  any  proposed  mechanism  must  fit.  In 
particular  the  attenuation  is  localized  in  the  upper  500  m  of  the  water 
column  and  it  must  be  very  insensitive  to  frequency.  Although  the  mag¬ 
nitude  of  the  assumed  attenuation  is  small,  it  does  not  differ  signifi¬ 
cantly  from  the  attenuation  in  sediments  at  these  low  frequencies.  This 
is  illustrated  in  Fig.  1 6 ,  where  the  attenuation  used  in  these  models  is 
compared  with  the  experimental  data  obtained  in  the  last  five  years  for 
oceanic  sediments.  These  attenuation  values  are  within  a  factor  of  ten 
of  each  other  at  50  Hz.  Since  the  path  length  in  water  is  often  ten 
times  greater  than  the  path  length  in  the  sediment  for  these  low  order 
modes,  the  attenuation  produced  by  the  waterborne  scatterers  has  the 
same  magnitude  as  the  attenuation  in  the  sediments. 

The  path  lengths  in  water  and  in  the  sediments  are  shown  in  Fig.  17 
as  a  function  of  bottom  grazing  angle.  These  curves  are  computed  for 
the  sound  velocity  profile  shown  in  Fig.  9  and  for  a  water  depth  of 
4883  m.  In  these  calculations  the  ocean  bottom  has  500  m  of  clay  over- 
lying  basalt.  The  velocity  contrast  between  the  water  and  clay  is  unity 
and  the  velocity  gradient  in  the  clay  is  1.23  sec"  .  The  path  lengths 
are  3een  to  differ  by  at  least  an  order  of  magnitude.  As  stated  above, 
these  differences  can  offset  the  differences  in  the  attenuation  of  the 
two  media,  and  thus  the  waterborne  scatterers  can  influence  the  mode 
attenuation  as  much  as  the  sediment  attenuation. 

In  summary  it  is  claimed  that  the  attenuation  observed  in  the  deep 
ocean  at  frequencies  below  100  Hz  is  produced  by  frequency  independent 
mechanisms  that  are  concentrated  in  the  upper  500  m  of  the  ocean. 
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